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In this review we will discuss the experimental and theoretical progresses in studying 
spin-orbit coupled degenerate atomic gases during the last two years. We shall first review 
a series of pioneering experiments in generating synthetic gauge potentials and spin-orbit 
coupling in atomic gases by engineering atom-light interaction. Realization of spin-orbit 
coupled quantum gases opens a new avenue in cold atom physics, and also brings out a 
lot of new physical problems. In particular, the interplay between spin-orbit coupling and 
inter-atomic interaction leads to many intriguing phenomena. Here, by reviewing recent 
theoretical studies of both interacting bosons and fermions with isotropic Rashba spin- 
orbit coupling, the key message delivered here is that spin-orbit coupling can enhance 
the interaction effects, and make the interaction effects much more dramatic even in the 
weakly interacting regime. 
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Many interesting phenomena in condensed matter physics occur when electrons 
are placed in an electric or magnetic field, or possess strong spin-orbit (SO) cou- 
pling. However, in the cold atom systems, neutral atoms neither possess gauge 
coupling to electromagnetic fields nor have SO coupling. Recently, by controlling 
atom-light interaction, one can generate an artificial external abelian or non-abelian 
gauge potential coupled to neutral atoms. The basic principle is based on the Berry 
phase effect 1^ and its non-abelian generalization^. An important application of this 
scheme is creating an effective SO coupling in degenerate atomic gases. Since 2009, 
Spielman's group in NIST has successfully implemented this principle and gener- 
ated both synthetic uniform gauge field ^, magnetic field electric field 1^ and SO 
coupling ^. We shall discuss the experimental progresses along this line in Sec. [l] 

The effects of SO coupling in electronic systems have been extensively discussed 
in condensed matter physics before and are also important topics nowadays. One of 
the most famous example is recently discovered topological insulators l^JiH In this 
review, we try to convey the message that SO coupling in degenerate atomic gases 
will bring out new physics which have not been considered before, mainly due to 
the interplay between SO coupling and the unique properties of atomic gases. For 
bosonic atoms, SO coupled interacting bosons is a system never explored in physics 
before. For fermionic atoms, since a lot of intriguing physics have been revealed 
during the last ten years by utilizing Feshbach resonance technique to achieve in- 
teraction as strong as Fermi energy, the interplay between resonance physics and 
SO coupling is definitely a subject of great interests. 
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A key point we want to emphasize in this review is that a nearly isotropic 
SO coupling will dramatically enhance the effects of inter-particle interactions, so 
that the interaction effects are not weak even in the regime where the interaction 
strength itself is small. This is because an isotropic Rashba SO couphng or a nearly 
isotropic SO couphng significantly changes the low-energy states of single particle 
Hamiltonian, as we shall discuss in Sec [2] In Sec. [3| we discuss many-body system 
of bosons. Because the single particle ground state has large degeneracy, it is the 
inter-particle interaction that selects out a unique many-body ground state and 
determines its low-energy fiuctuations. In Sec. |4j we discuss many-body system of 
fermions. Because the low-energy density-of-state (DOS) is largely enhanced, the 
interaction effects become much more profound, in particular for weak attractions. 
A brief summary and future perspective are given in Sec. [5] 

1. Synthetic Gauge Potentials and Spin-Orbit Coupling 

In 2009, Spielman's group in NIST first realized a uniform vector potential in Bose- 
Einstein condensate (BEG) of ^^Rb ^. In this experiment, two counter propagating 
Raman laser beams couple |F, mp) = ll^ — 1) level of ^^Rb to |1, 0) level, and couple 
|1,0) level to |1, 1) level, as shown in Fig. [TJa) and (b), which can be described by 
the Hamiltonian 



where ko = 27r/A, A is the wave length of two lasers. 2ko is therefore the momentum 
transfer during the two-photon processes, ci = Ai+(5a; + A2, and 62 = Ai+(5cj — A2, 
where Ai denotes the linear Zeeman energy, Su denotes the frequency difference of 
two Raman lasers, and A2 is the quadratic Zeeman energy. 

Applying a unitary transformation to wave function <l> = /7^, where 



When both ei and 62 are large, the single particle energy dispersion of i^eff is shown 
as Fig.jljc), which displays a single energy minimum at finite kx- In this regime the 
low energy physics is dominated by a single dressed state described by ^{kx—Ax)'^^ 
where Ax is a constant. This leads to a uniform vector gauge field. 

In a follow up experiment, Spielman's group applied a Zeeman field gradient 
along y direction to this system ^ . In this case. Ax becomes a function of y instead 




(1) 




and the effective Hamiltonian becomes 




(3) 
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Fig. 1. (a) A schematic of NIST experiment, in which two counter propagating Raman beams are 
apphed to ^^Rb BEC. (b) A schematic of how three F = 1 levels are coupled by Raman beams, (c) 
Dispersion in the regime of uniform vector potential, (d) Dispersion in the regime of non-abelian 
gauge field. 

of a constant. It gives rise to a non-zero synthetic magnetic field Bz = —dyA^ 7^ 0. 
They observed a critical magnetic field above which many vortices are generated in 
the BEC ^ . In another experiment ^, they made time dependent which gives rise 
to a non-zero electric filed = —dtA^ 7^ 0. They observed collective oscillation of 
BEC after a pulse of electric field ^. 

By tuning the Zeeman energy and the laser frequency, one can also reach the 
regime where Ai + (^cj ~ A2, and thus 62 ~ 0, while ei ~ 2A2 is still large. In 2011, 
Spielman's group first reached this regime and showed that a SO coupling can be 
generated^. As shown in Fig.jljd), in this regime the low-energy physics contains 
two energy minima which are dominated by |1,— 1) and |1,0) states, respectively. 
Therefore we can deduce the low-energy effective Hamiltonian by keeping both 
|1, — 1) and |1,0), and rewrite the Hamiltonian as 



where h = €2. Similarly, by applying a unitary transformation to the wave function 
with 
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one reaches an effective Hamiltonian that describes SO coupling 
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In fact, upon a pseudo-spin rotation —cfz and cfx^ the above Hamiltonian 

is equivalent to 

where the first term can be viewed as an equal weight of Rashba (kx^x + ^y^y) 
and Dresselhaus (kx'^x — ^y^y) SO coupling |^ The Hamiltonian of Eq. [7|can also 
be viewed as a Hamiltonian with synthetic non-abelian gauge field, since the vector 
potential = — /cqct^ does not commute with the scale potential ^ = ^cix -\- ^cr^. 

Later on, Campbell et al. discussed how to generalize the NIST scheme to create 
SO coupling in both x and y direction, and nearly isotropic Rashba SO coupling tl^. 
Xu and You recently introduce a dynamics generalization of the NIST scheme 
Sau et al. discussed an explicit method to create Rashba SO coupling in fermionic 
^^K system at finite magnetic field, where a magnetic Feshbach resonance is avail- 
able Beside the NIST scheme, there are also other theoretical proposals that use 
A- type or tripod system to generate synthetic magnetic field EIEIEIIH^ non-abelian 
gauge field with a monopole and SO coupling 'H*i222*£ll However, for those pro- 
posals using dark states IMlllIElIMllI, collisional stability is a concern since there 
is always one eigen-state whose energy is below the dark state manifold, and multi- 
particle collision can lead to decay out of the dark state manifold. A recent review 
paper by Dalibard et al. has discussed different proposals in detail 



2. Single Particle Properties with Rashba Spin-Orbit Coupling 

In the rest part of this review we will consider SO coupling in both x and y directions, 
whose Hamiltonian is given by 

i^o = ^ [{kx - nxCTxf + {ky - nyCFyf + kl] (8) 

and in particular, we will consider the most symmetric Rashba case — — 
n > where the physics is most interesting. In this case, the Hamiltonian can be 
rewritten as 

H^ = ^ (ki - 2n\^^ -G^n^^kl) (9) 

Obviously, spin is no longer a good quantum number. However "helicity" is a good 
quantum number. "Helicity" ± means that the spin direction is either parallel or 
anti-parallel to the in-plane momentum direction. For these two helicity branches, 
their dispersion are given by 

ek± = ^{kl^ 2Kk^ + + kl) (10) 



^In many literature, Rashba SO coupling denotes kxCTy — kyax, while Dresselhaus SO coupling 
denotes kxdy + kydx- They are equivalent to the notations used in this paper by a pseudo-spin 
rotation ax — ^ ~'^y cind ay ax- 
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Fig. 2. Energy dispersion with kz = (a) and density-of-state (b) for Rashba spin-orbit coupled 
single particle Hamiltonian. In (a), "helicity" is " + " for red solid line and is " — " for blue dashed 
line. E'r = ^c^/(2m) is introduced as energy unit. 



where = (kxjky) and k± = yk'^-\-ky. This Hamiltonian displays a symmetry 
of simultaneously rotation of spin and momentum along z direction. 

As shown in Fig. |2|a), helicity plus branch has lower energy. The single par- 
ticle energy minimum has finite in-plane momentum k± = k,^ and all the single 
particle states with same k± = and kz = but different azimuthal angle are 
degenerate ground states. The single particle DOS also has non-trival feature. In a 
three-dimensional system without SO coupling, DOS vanishes as ^/e at low ener- 
gies. However, in this system, as shown in Fig. |2|b), the low-energy DOS becomes 
a constant when e < £^r, similar to a conventional two-dimensional system. 

Without SO coupling, the single particle Hamiltonian has a unique ground state 
at k = 0. At zero-temperature, bosons are all condensed at k = state. If the inter- 
action is weak, interaction effect is perturbative which only creates finite quantum 
depletion to the zero-momentum condensate. Also, because the vanishing DOS at 
low energies, two-body bound state appears only when the attractive interaction is 
above a threshold. Far below the threshold, when the attractive interaction is weak, 
the strength of fermion pairing and the fermion superfiuid transition temperature 
are both exponentially small. That is to say, without SO coupling, not surprisingly, 
the interaction effect is weak in the regime where the interaction strength is small. 

With SO coupling, the effects of interaction are significantly enhanced even when 
its strength is still weak, because SO coupling significantly changes the single par- 
ticle behaviors as discussed above. First, because the single particle ground state 
is not unique now, the ground state of a boson condensate is also not unique if 
there is no interactions. That is to say, it is the interaction that selects out a unique 
ground state among many possibilities. In this sense, the effect of interaction is 
non-perturbative even for very weak interactions. Secondly, because the DOS is a 
constant at low energies, a two-body bound state appears for any weak attractions. 
Therefore, for weak attractive interactions, both pairing gap and superfiuid transi- 
tion temperature are largely enhanced by SO coupling. In particular, the superfiuid 
transition temperature can reach the same order as the Fermi temperature even for 
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weak attractions. These points will be discussed in the next two sections in more 
detail. 

On the other hand, we shall always keep in mind that in real experiment, one can 
never achieve a perfect isotropic Rashba SO coupling. Therefore, in our following 
discussions, we always first obtain interesting results in the isotropic Rashba limit, 
and then discuss how robust the results are if there is a little anisotropy, namely, a 
slight difference between tZx and hZy. 



3. Spin-Orbit Coupled Bose Gases 

We will first consider (pseudo-) spin- 1/2 bosons. We shall first consider the mean- 
field ground state and then discuss the effects of quantum and thermal fluctuations 
on top of mean- field saddle points. Let us first consider the most simplified form of 
interactions 

^int = J d^r{gnl{r)^gnl{r)^2gi2ni{r)n2{r)) (11) 

By minimizing Gross-Pitaevskii energy functional, Wang et at. found that the mean- 
field ground state has two different phases depending on the sign of gi2—g^^- When 
912 < the system is in the "plane wave phase" , where all bosons are condensed into 
a single plane wave state, and the direction of plane wave is spontaneously chosen 
in the xy plane. For instance, if the plane wave momentum is along x direction, the 
condensate wave function is given by 



^=A/^e-M . (12) 




The density of each component is uniform, but their phases modulate from zero 
to 27r periodically, as found from numerical solution of Gross-Pitaevskii equation 
and shown in Fig. [sj^a). This state spontaneously breaks time-reversal, rotational 
symmetry and the U (1) symmetry of superfluid phase. When gi2 > g^ all bosons are 
condensed into a superposition of two plane wave states with opposite moment ums, 
whose condensate wave function is given by 



In this phase, the spin density ni — n2 = pcos{2hzx) which has a periodic modulation 
in space, as shown in Fig. [sjb), and therefore is named as "stripe superfluid". The 
direction of the stripe is also spontaneously chosen in the xy plane. Here, without 
loss of generality, we choose it along x direction. In this state, the high density regime 
of one component coincides with the low-density regime of the other component, so 
that the inter-component repulsive interaction is maximumly avoided. That is the 



^ There are also proposals of non-mean-field fragmented state, like NOON state, as ground state of 
SO coupled bosons however, the conventional wisdom is that such a state is very fragile when 
external perturbations are present. 
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Fig. 3. (a) Phase of condensate wave function in the "plane wave phase"; (b) Spin density in the 
"stripe superfluid" phase. 

reason why the spin stripe state is favored when gi2 is larger than g. In addition to 
the U{1) superfluid phase, this state also breaks the rotational symmetry (but keeps 
the reflection symmetry), and the translational symmetry along the stripe direction. 
Hence, symmetry wise, it can also been called "smectic superfluid". In contrast to 
the "plane wave phase" , this state does not break time-reversal symmetry. 

In principle, the condensate wave function can be any superposition of single 
particle ground states as 



One may wonder why only a single state or a superposition of a pair of states 
are favored by interactions. In fact, if there are more than a pair of single particle 
states in the superposition, the condensate wave function will exhibit interesting 
structure, such as various types of skyrmion lattices. And if all the degenerate 
states enter the condensate wave function with equal weight and specified relative 
phases, condensate will exhibit interesting structure of half vortices, as first proposed 
by Stanescu, Anderson, Galitski 1^ and Wu, Mondragon-Shem However, such a 
state is not energetically favorable in spin- 1/2 case for a uniform (or nearly uniform) 
system. This is because that the interaction part can be rewritten as 



where n(r) = ni(r) + n2(r) and Sz = ni(r) — n2(r). The 5'^-term can be satisfied 
by either choosing the "plane wave phase" or the "stripe phase", while the n^- 
term with positive coefficient always favors a uniform density. One can easily show 
that, if there are more than a pair of states in the superposition, the condensate 
density will always have non-uniform modulation, which causes energy of n^-term. 
Similar situation has also been found for spin-1 Hamiltonian However, there 

are several situations where skyrmion lattices or half vortices are found as ground 
state. When a strong harmonic confinement potential V^(r) = mijor\/2 is applied 
to the system, condensate density is no longer uniform because of the trapping 
potential and the requirement from n^-term becomes less restrictive. In addition. 
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(15) 
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one notes in the limit of zero interaction, the ground state in a harmonic trap can 
be solved exactly and it is a half- vortex state EHHEH Recently, three groups have 
found that if one largely increases the confinement potential so that a = ^/h/muJ is 
comparable to 1//^, or reduces the interaction energy to be comparable to hu^ the 
ground state will evolve continuously to skyrmion lattice phases, and finally to half 
vortex phases mmEH Besides, Xu et al. and Kawakami et al. found in spin-2 case, 
because of an addition interaction term (which favors cyclic phase in absence of SO 
coupling P ^l^^ l ) , there exists a parameter regime where various types of skyrmion 
lattice phases are ground state even for a uniform system. 

Back to the discussion of a nearly uniform spin- 1/2 systenj^ the "plane wave 
phase" and the "stripe phase" are in fact two very robust mean-field states. In 
real situation, the interactions between two pseudo-spin states have much more 



complicated form than the simplified form of Eq. (11). Both Yip and Zhang et al. 
considered a specific type of complicated interaction form using a concrete real- 
ization of Rashba SO coupling, and they found the ground state is still either the 
"plane wave phase" or the "stripe phase" And also, the SO coupling is always 

not perfectly isotropic, say, Kx > t^y^ then the Hamiltonian itself does not have rota- 
tional symmetry anymore. The single particle energy has two minima at 0, 0) 
instead of a continuous circle of degenerate states. At mean-field level, the effect of 
anisotropic SO coupling is to pin the direction of plane wave or stripe into certain 
direction {x direction for this case). The NIST experimental situation discussed in 
Sec. [T] corresponds to the case K.y = and also with a Zeeman field. As shown by 
Ho and Zhang and also in the experimental paper the phase diagram of this 
system also only contains such two phases. 

To go beyond the mean-field description, three different approaches have been 
tried so far. The first is the effective field theory approach which treats Gaussian 
fiuctuations of low-energy modes. The second is Bogoliubov approachEHmm which 
more focuses on the gapless phonon excitations. And the third is the renormalization 
approach ESE2II2I, which discusses how scattering vertices are renormalized by high 
order processes. These different approaches address beyond-mean-field effects from 
different perspectives, and the results are consistent with each other where they 
overlap. 

Taking effective field theory approach as an example, for the "stripe" phase, the 
superfiuid phase and the relative phase u between two momentum components 
are two low- lying modes: 

In fact, the relative phase u describes the phonon mode of stripe order. The dy- 



^For typical experimental parameters, the BEC is in the nearly uniform regime rather than strong 
harmonic confinement regime. 
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namics of and u fields are governed by an effective energy function!^ 



njST _ P 



u) 



(17) 



where a > 1 is a constant. For the "plane wave" phase, the super fuid phase is the 
only low energy mode 



^pw = v^eH-^"^ y ^ J (18) 
and its effective energy is derived a^^ 



(19) 



One notes that in the "stripe" phase, the quadratic term {dyuY is absent in Eq. 
and in the "plane wave" phase, the quadratic term {dyOY is absent in Eq. 



. This is in fact a manifestation of rotational symmetry in this system. Similar 
effective theory has also been found for "FFLO" state in fermion superfluid ElEH 
Such an energy function is also the classical energy of smectic liquid crystal. 

In two dimensions, the finite temperature phase transition is driven by prolif- 
eration of topological defects. For usual XY model, both the energy of topological 
vortex and its entropy logarithmically depend on system size. Thus, only above 
a critical temperature, entropy wins energy and the topological defects proliferate 
which drives system into normal phase. However, in the "stripe" phase, because the 
absence of {dyuY term, the energy for a topological defect of li, i.e. a dislocation in 
the stripe, no longer logarithmically depends on system size. Hence it will lose to en- 
tropy at any finite temperature, and the proliferation of these dislocations will melt 
the stripe order. This restores translational symmetry and drives the system into a 
boson paired superfluid phase . Such a boson pairing state can also be predicted 
by looking at pairing instability of normal state with renormalized interactions 
By considering the renormalization effects of scattering vertex from high energy 
states, Gopalakrishnan et al. ^ and Ozawa, Baym ^^9^ found that the scattering 
amplitudes between two states with opposite momentum become smaller and even 
vanishing, which makes the "stripe phase" more favorable at the low-density limit, 
and meanwhile leads to more significant the fluctuation of stripe order For same 
reason, in the "plane wave" phase, the superfluid phase will immediately disorder 
at finite temperature and the system becomes normal. 



In addition, the effective theory Eq.(17) and Eq. (19) also imply that there is a 



Goldstone mode which has linear dispersion along x direction (direction of stripe 
or plane wave momentum), and quadratic dispersion along y direction (direction 
perpendicular to the direction of stripe or plane wave momentum). Same results 
have been reached by Bogoliubov calculation ^ElMl Similar behaviors of Goldstone 
modes also exist in "FFLO" phase of fermion superfluid E1EII3H 
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When the SO couphng is anisotropic {tvx there is no rotational symmetry 

in the Hamiltonian. Therefore, one wih have (pyU)^ in the stripe phase, and {dyOY 
in the plane wave phase. However, the coefficients of those terms are propositional 
to 1 — (tZy/hZx)'^' In the regime hZy/hZx is very close to unity, the energy of a disloca- 
tion is much smaller than the energy of a vortex or a half vortex. Hence, the system 
will undergo two Kosterlitz-Thouless phase transitions. At a lower critical tempera- 
ture, dislocations proliferate and the system becomes paired superfluid. Then, at a 
higher critical temperature, vortices proliferate and the system becomes normal. A 
completed phase diagram in term of interaction parameter, SO coupling anisotropy 
and temperature is given by Jian and ZhailSIl. 

For SO coupled bosons, many questions remain open. Recently several works 
begin to address the questions about the superfluid critical velocity vortices in 
presence of external rotation EIEIEII, the effects of dipolar interaction^^, the collec- 
tive modes superfluid to Mott insulator transition in a lattice the interplay 
between magnetic field and SO coupling and the dynamical effects nearby Dirac 
point due to SO coupling | 54 | 55 | 56 [ r^^ie research along this direction will definitely 
reveal more interesting physics and stimulate more interesting experiments. 



4. Spin- Orbit Coupled Fermi Gases across a Feshbach Resonance 

Even for a non-interacting system, the thermodynamic behavior of a Fermi gas is 
dramatically changed by a strong SO coupling because of the change of the low- 
energy DOS . In this section, we focus on Fermi gas with attractive interaction, 
and in particular, across a Feshbach resonance. The interaction part is modeled as 



Hi, 

where 



= gl dW^{r)i,l{r)M^)Mr) (20) 



1 m ^ 1 . 



Here we relate the bare interaction g to 5- wave scattering length via Eq. (21 ), which 
is the same as the scheme widely used in free space without SO coupling. This is 
based on the assumption that the interaction Hamiltonian is not changed by SO 
coupling, and the here should be understood as scattering length in free space. 

Using this interaction Hamiltonian, Vyasanakere and Shenoy ffist studied the 
two-body problem across a Feshbach resonance Because the low-energy DOS is 
now a constant, an arbitrary weak attractive interaction will give rise to a bound 
state. Similar situations are two-body problem in two-dimension and Cooper prob- 
lem in three-dimension in absence of SO coupling, where DOS are also constants. 



^This is equivalent to assume that the s-wave pseudo-potential is still a valid approximation 
for a short-range realistic potential in presence of SO coupling. In fact, the validity of such an 
approximation is not quite obvious, and it has only been examined recently by Cui*^. 
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Fig. 4. (a) Gap size A as a function of ^t/Zcp for different values of l/(/cFas). (b) Size of Cooper 
pair in xy plane and along z direction Iz as a function of ^c/Zcp at resonance with as = oo; (c) 
Superfluid transition temperature Tc/Tp as a function of l/(/cFas) for different values of hi/k-p. 
Reprinted from arXiv: 1105.2250 (Phys. Rev. Lett, to be published). 

The binding energy can be easily calculated by looking at poles of T-matrix 
or by reducing the two-body Schrodinger equation to a self-consistent equation 
The two-body properties at the BCS side and at resonance regime are significantly 
changed by SO coupling. At weakly interacting BCS side with small negative tza^^ 
the binding energy behaves as 

h^K^ 4 

E^^-- -e -i-i. 22 

2m 

where a large binding energy can always been reached by increasing the strength 
of SO coupling. At resonance with = oo, is the only length scale in the 
two-body problem and therefore one has a universal result 

^b = -0.88^. (23) 

While for the BEG side with small positive A^ag, at leading order E]j is still given by 
/{2mal) and is not affected by SO coupling. Moreover, because of SO coupling, 
the two-body wave function has both singlet and triplet components. For a two- 
body bound state with zero center-of-mass momentum, the wave function behaves 

^ = Vs(r)i n - ;t> + ^:(r)i n) + u) (24) 

where V^sCr) and V^alr) are symmetric and anti-symmetric functions, respectively. 
Furthermore, by looking at binding energy with finite center-of-mass momentum, 
one can determine the effective mass of molecules (two-body bound state). Hu et al. 
1^ and Yu and Zhai found that at the BCS limit, the effective mass of molecule 
finally saturates to 4m, and at resonance, the effective mass is a universal number 
of 2.40m. In the BEC limit, the effective mass saturates to 2m as conventional case. 

Generalizing conventional BEC-BGS crossover mean-field theory to the case 
with SO coupling and equal population (n^) = (n^), one can show that the system 
remains gapped for all /cf<^s7 although there are triplet p-wave components 
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and the pair wave function obtained from the mean-field theory finahy approaches 
the wave function of two-body bound state (i.e. molecule wave functionjiniH Hence, 
it is still a crossover as l/(/cFas) changes from negative to positive. However, on the 
other hand, the change of low-energy DOS and the presence of two-body bound 
state at the BCS side and resonance regime will significantly change the properties 
of crossover |6Q|6i|62(64[ example, as shown in Fig. |4|a), the pairing gap at the 
BCS side is dramatically enhanced when K/k^ is comparable or larger than unity j^] 
For such a strong SO coupling, the DOS at Fermi energy becomes a constant, and 
is much larger than the DOS in absence of SO coupling. That is the reason why 
the pairing effects become much dramatic even for same interaction strength. For 
another example, in absence of SO coupling, Fermi energy is the only energy scale at 
resonance, and therefore the size of Cooper pair is a universal constant times l/k-p- 
SO coupling introduces another scale at resonance, which is For large tz^ as 
the pairing gap approaches two-body binding energy, the size of Cooper pairs also 
approaches I/k. Fig. |4]^b) shows that the behavior of / undergoes a crossover from 
l/kp to as f<i/kF increases. This plot also shows that the size of Cooper pair in 
the xy plane is different from the size along z direction, namely, the Cooper pairs 
are anisotropic. In addition, one can also show that the superfiuid transition 

temperature at the BCS side can also be enhanced a lot by SO coupling. For large 
enough SO coupling, it eventually approaches the BEC temperature of molecules 
with mass 4m EaMI, which is a sizable fraction of Fermi temperature. The critical 
temperature across resonance is first estimated by Yu and Zhai as shown in Fig. 

Be). 

If SO coupling is slightly anisotropic, DOS at very low-energy will finally vanish. 
However, there is still a large energy window < h^H? /{2m) where the DOS is greatly 
enhanced by SO coupling. Hence, pairing gap will still be enhanced as long as the 
density of fermions is not extremely low. Moreover, although it is no longer true 
that arbitrary small attraction can cause a bound state, the critical value for the 
appearance of a bound state will move from unitary point ttg = co to the BCS side 
with negative Once the bound state is present, it will infiuence the universal 

behavior of pair size at resonance and superfiuid critical temperature as discussed 
above. 

For the imbalanced case, the phase diagram becomes more richer. Several groups 
have studied the phase diagram in presence of a Zeeman field |65|66|67|68{69] 
in various other circumstance l^"r^r^r^r^r^r^r^r^r^l^"l[l They have shown that. 



^ In cold atom system, SO coupling is generated by atom-light coupling, and therefore n is on the 
order of the inverse of the laser wave length. And since the laser wave length and the inter-particle 
distance are comparable (between ~ 0.1/im and ~ Ijim) in atomic gases, the strength of SO 
coupling in cold atom systems can naturally reach the regime ^t/Zcp ~ 1. 

^ Iskin and Subasi studied SO coupled Fermi gas with mass imbalance They use mixture of 
different species as motivation of this study. We caution that the current way of generating SO 
coupling is based on light coupling of different internal state of atoms, which can not generate SO 
coupling if different internal states are different atomic species 
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instead of a crossover, there are phase transitions between topological and non- 
topological phases | 65 | 66 | 67 | 68 [ xhey have also discussed how SO coupling influences 
the competition between a uniform superfluid and a phase separation [MIIZI 

5. Summary and Future Developments 

SO coupled quantum gases with interactions are new systems in cold atom physics. 
Moreover, SO coupled bosonic system has never been thought in physics before. 
Currently our understanding of this system is still very limited, and many questions 
remain open. However, even from our limited experience with this new system, one 
can already get a feeling that this system has many unusual behaviors. This gives 
a lot of opportunities for theorists and experimentalists. 
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